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SUMMARY 


A method is presented for determining the control gearing and time 
lag necessary for a specified damping of the motions of an aircraft 
equipped with an autopilot having constant-time-lag characteristics. 

The method is applied to a typical present-day airplane equipped with an 
autopilot which applies rudder control proportional to the yawing angular 
acceleration. The types of motion predicted for this airplane -autopilot 
system by this method are In very good agreement with the airplane 
motions calculated by a step-by-step procedure. 

For some values of control gearing there may exist more than one 
range of time, lag for which the motions of the assumed aircraft- 
autopilot system will have more than a specified amount of damping. 

INTRODUCTION 


Many present-day aircraft designed to fly at transonic and super- 
sonic speeds have exhibited poor lateral stability characteristics. As 
a result, much interest has been shown in automatic stabilization systems 
as a means of obtaining satisfactory stability for high-speed flight. 

In analyzing the effect of a particular autopilot, the usual practice 
has been to determine whether the aircraft -autopilot system is stable 
by employing methods such as those presented in reference 1 (frequency- 
response analysis) and reference 2 (Nyquist plots) . For autopilots 
characterized by constant time lag ( linear variation of control gearing 
and phase lag with frequency) the same information can be obtained by 
carrying out an analysis such as that of reference 3* This type of 
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information is important in the analysis of an autopilot-aircraft 
system hut affords little quantitative indication of the degree 
of stability which can be obtained by use of a particular autopilot. 

Since one of the purposes of equipping an aircraft with an autopilot is 
to improve its stability, the determination of the type of autopilot 
frequency response that will result in satisfactory stability appears to 
be of greater importance. The degree of stability which can be obtained 
by equipping an aircraft with a linear autopilot can be calculated by 
the methods of reference 4, but, if an over-all picture is desired of the 
variation of stability for different combinations of control gearing and 
constant time lag, the calculations become too laborious for practical 
application. A method of obtaining this information was discussed very 
briefly in the appendix of reference 4, but has been found inadequate 
for a comprehensive analysis of autopilot systems which exhibit constant- 
control-gearing and constant -time -lag characteristics. 

The purpose of the present paper, therefore, is to extend the 
concept discussed in the appendix of reference 4 and thereby present a 
rigorous method of obtaining the combinations of constant control gearing 
and constant time lag necessary to provide specified amounts of damping , 
up- to the maximum obtainable, to the motions of a linear oscillating 
system, with application to the automatic stabilization of aircraft. 

The method presented is applicable, in a strict sense, only to control 
systems which exhibit linear frequency-response characteristics, but, 
under certain conditions discussed in the paper, the analysis is valid 
for systems having a frequency response similar to that of a constant- 
time-lag system over a limited range of frequencies. 

Although in this paper the method is applied to the problem of 
automatic stabilization of aircraft, it is, in general, applicable to 
any linear oscillating system with constant time lag. 


SYMBOLS AND COEFFICIENTS 

0 angle of roll, radians 

angle of yaw, radians 

P angle of sideslip, radians (v/V) 

r yawing angular velocity, radians per second (di/dt) 

V yawing angular acceleration, radians per second per second 

( d 2 ^/dt 2 ) 
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k Xo 

k Zo 

KXo 

K Zo 

K X 

K Z 


rolling angular velocity, radians per second (d0/dt) 
sideslip velocity along lateral axis, feet per second 
airspeed, feet per second 
mass density of air, slugs per cubic foot 

ipV^) 

wing span, feet 

wing area, square feet 

weight of airplane, pounds 

mass of airplane, slugs (w/g); any integer 
acceleration due to gravity, feet per second per second 
relative-density factor (m/pSb) 

inclination of principal longitudinal axis of airplane with 
respect to flight path, positive when principal axis is 
above flight path at nose, degrees 

angle of flight path to horizontal axis, positive in a climb, 
degrees 

radius of gyration in roll about principal longitudinal axis, 
feet 

radius of gyration in yaw about principal vertical axis, feet 

nondimensional radius of gyration in roll about principal 
longitudinal axis ^kx Q jb^ 

nondimensional radius of gyration in yaw about principal 
vertical axis ^k^b^ 

nondimensional radius of gyration in roll about longitudinal 
stability axis ^\/ Kx o 2cos2 ^ + K 2 c 2sin ^ r l 

nondimensional radius of gyration in yaw about vertical 
stability axis (^K^^cos^rj + K^^sin^T] 




dynamic pressure, pounds per square foot ^ 
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K xz 

nondimensional pro duct -of - 

inertia parameter 


((kZq 2 - KXo 2 )sin r, cos 

r ‘) 

c L 

/w 

trim lift coefficient f — 

cos y\ 
qS ) 

c z 

rolling -moment coefficient 

/Rolling moment \ 
V qSb ) 

c n 

yawing -moment coefficient 

/Yawing moment \ 
\ qSb ) 

c Y 

lateral-force coefficient 

/Lateral force \ 

\ qS y 


Cz P ■ w 



dC n 

33 



3Cy 

W 


3c 


n 


^n T 


32^- 

2V 



j n p 


= 3Cn 

3^ 

°2V 


C, = ^ 
P 32^ 
°2V 


C Y =^I 

P 32^ 

°2V 


C Y. = 


c3Cy 

^rlD 

% 
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Sb„ 


'Qb^ = 


. = ^c n 

" n 5a ' 5b^ 


C Z5 = |£i 
6a Sb a 


^5 r 

%> 

C Y8 r 

t 

Sb 

Db 

P 

T l/2 


8 r> 8 a 


a 

co 

(Dg 


Sc 

Sb^ 

Sc^ 

Sbl 


z 


SCy 

Sbl 


time, seconds 

nondimensional time parameter based on span (Vt/b) 
differential operator /— i 

V ds b/ 

period of oscillation, seconds 

time for amplitude of oscillation to damp to one-^half its 
original value, seconds 

deflection of control surfaces, radians 

real part of complex root of characteristic stability equation 
angular frequency, radians per second 
nondimensional angular frequency 

complex root of characteristic stability equation (a + io> s ) 
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time lag between signal for control and its actual motion, 
seconds 


nondimensional time 


lag (b) 


amplitude of control-surface oscillation produced by autopilot 
in response to oscillation of airplane acceleration 

/IS, 


V 


autopilot 


) 


k s = 


K 2 t 


= ¥§k 
b 2 


phase angle, radians 


EQUATIONS OF MOTION 


The nondimensional linearized equations of motion, referred to the 
stability axes, are: 

Yawing 

•2|»b( K zV* + %Z »&) * V + 1% M + ¥nr “b* + C n 6r 6 r (la) 

Rolling 

2pb( K X 2l} b 2 0 + K XZ Db 2 ^) = c ZpP + ^ + i°ir ^ + Cz 6 a &a 

Sideslipping 

2n b (DbP + Db*) = C Yp p + |c Yp Db0 + C L 0 + |c Yr Db* + (C L tan 7 )* Uc) 

The derivatives Cng a > CYg a > an( t CY§ r have been neglected 

in equations (l). 

If a constant-time-lag autopilot which applies rudder control pro- 
portional to the nth derivative of the yawing displacement is installed 
in the airplane, the equation for 5 r as a function of s b becomes 
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^r( s b) 


d& r 

^Db n ^ 


I>b n t(s t - t s ) 


( 2 ) 


where the term Dt. i(sh - t s ) signifies that the control deflection at 

a given instant is proportional the the nth derivative of which 
e ^ sted at a fixed time t s previous to the given instant. The term 

is the c °ntrol-gearing ratio of the autopilot. Similarly, if a 


constant-time-lag autopilot which applies aileron control proportional 
to the nth derivative of the rolling displacement is installed in the 
airplane, the equation for 6 a as a function of s b becomes 


Sa<Sb) ‘s5^ D> ’ ,10(sb " Ts) 


( 3 ) 


For purposes of illustration, assume that the type of autopilot 
described by equation (2). is to be investigated. Therefore, the value 
of 5 r given by equation (2) is substituted for & r in equations (l). 
The aileron deflection 6 a is assumed zero. When 0 o e^ S ^ is substi- 
tuted for 0, \J/ 0 e*' & *- ) for i)r, and p Q e^ Sl3 for (3 in the resultant 
equations written in determinant form, A must be a root of the charac- 
teristic stability equation 

AX. 5 + baA + CA 3 + DX 2 + EX + k s X n e~ Ts \A' A 3 + B 1 A 2 + C' A + D' ) = 0 (4) 

where A, B, C, D, E, A f , B* , C f , and D* are functions of the mass and 
aerodynamic parameters of the airplane. The expressions for A, B, C, 

D, and E are given in reference and, for this particular case, 

A* = -^ b 2 K X 2 Cn 8r 

B' = ^ 2 ^ b K^ 2 CYp + hb 2 Zp)^n5 r 

+ £% C *p) Cn &r 

D ' = C ^°n 5r • 

Ss r 

and k q = . 

a Db n i 
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The time to damp and the period of the lateral oscillation in seconds 
are, given "by the expressions 


*1/2 


p — b 
CDo V 


( 5 ) 


where a and <u s are the real and imaginary parts, respectively, of a 
complex root of equation ( 4) . 


ANALYSIS 

The following assumptions are made in the derivation of the method 
which may he used to determine the combinations of k s and t s neces- 
sary to provide a specified amount of damping of the lateral oscillations 
of an automatically stabilized airplane: 

(l) The characteristic stability equation of an airplane equipped 
with an autopilot which applies rudder control proportional to the nth 
derivative of the yawing displacement which existed at time s^ - t s is 
given by equation (*0; that is, the autopilot is assumed to have linear 
characteristics. 


(2) This stability equation has as one of its roots X. = a + i<n s . 
(This root X represents an oscillatory mode of motion, the time to 
damp to one-half amplitude and the period of which are given by equa- 
tions (5) • ) 

The stability equation of the airplane-autopilot system (equa- 
tion (4)) may be rewritten in the following form: 

- T S X. _ A\5 + BX 1 * + CX 3 + DX 2 + EX , 

SS -X*(a'X 3 + B’ X 2 + C'X + D’) 


If X = a + iois is substituted in each side of equation (6), the 
condition that X be a root of the characteristic’ stability equation is 
that the complex number Ai + iB]_ obtained from the left-hand side must 

be equal to the one obtained from the right-hand side A 2 + iB2- The - 
quantities. Ai + iBx and A 2 +. iB 2 may be represented by the expres- 
sions Eie iQl and R 2 e i02 , respectively. Therefore, this requirement 
is equivalent to saying that R]^ = R 2 and 0]. = 6 2 if X is to be a 
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root of equation (4). If X = a + ico s is substituted in the left-hand 
side of equation (6), the following expression results: 

k s e " Tsa /cos T s co g - i sin TgOiA = + iB 1 = R-^e^h 


Since the complex number A]_ + iBj is unchanged when its phase angle is 
increased by integral multiples of 2it,. the following identities exist: 


E i s (*.) « 

m 


■(t») 


in 


6]_ s 2jtm - ( t g ) <Ug 


where m - 0, 1, 2, . . . . If X = a + icn s is substituted in the 

right-hand side of equation (6), the values of R 2 and ' 6^ can be 

calculated numerically. For each value of m, the value of (t,,) 

a m 

necessary to make 8 ^ equal to- 8^ can be calculated and since a is 

a known fixed quantity, the values of (k ) necessary to make R, 

B m -L 


equal to R 2 
and (k B ) 

X = a + ioDg 


can also be determined. For the combinations (r a ) 

° 1 

(t s )o (^s)p^ • # • ( T s ) and (k 3 ) , therefore, 

m m 

is a root of the characteristic stability equation (equa- 


tion (4)). The parameter m in the preceding equations is necessary 
if the analysis is to hold for values of a which represent dampings 
less than the natural airplane damping . More specifically, the ranges 
of t s and k s which result in greater stability than the natural 

stability of the airplane alone cannot be obtained correctly by 
utilizing the methods as presented m reference 4. For example, the 
curves presented in figure 10 of reference 4 fail to provide the 
important information that, for small values of j the autopilot- 


airplane system discussed in that reference is less stable than the 
airplane alone. In addition, the parameter m is necessary if the 


analysis is to hold for values of These values of T must 

«>s 3 

be included in a theoretical analysis of the problem since, as t 8 is 

p -rr 

increased by multiples of — > the phase relation again exists which 

' 8 

is necessary for the system to have an oscillatory mode the, damping of 
which is indicated by a. The analysis of reference 4 considered only 
the case of m = 1 and was therefore restricted to values of r a <^r- 

S CSg 

or,, in effect, to values less than the natural period N of the oscillation. 
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For a given value of a corresponding to some desired value of 
T^/2 0:f the lateral oscillation, the analysis may be made throughout 

the range of co s and the required combinations of (k s ) m and ( T s ) m 
determined for each o> s . For each m, therefore, a curve representing 
this value of a may then be plotted as a function of (k s ) m and (t s ) 

This family of m curves; will divide- the k s ,T s plane into regions for 
which the lateral oscillations will, have more or less damping than that 
indicated by the value of a being investigated. 

General expressions for (T s ) m and (k s ) m for each oo s , in terms 
of R 2 and 62 , are as follows: 


(T S ) m = (r s ) n + 


2jnn 
ts 0 


(ks)* = (ks) 0 e 


’s 

2jona/cn s 


where 


< T =>o “-fir 


(ks)g - Rp e 


( T s ) 


0 a _ 


= R 0 e 


-(^A^s 


Thus, for any value of . o) S , it is necessary only to obtain (t s )q and 
(k s )g since, by use of the preceding generating functions, the required 


combination of 
of m. 


t s and k s are readily obtained for any integral value 


ILLUSTRATIVE EXAMPLES AND DISCUSSION 


Since the purpose of installing an autopilot is to improve the 
lateral stability of the airplane, one logical value of a to investi- 
gate would be that value which corresponds to the damping of the lateral 
oscillation of the airplane' without the autopilot installed, hereinafter 
referred to as airplane damping. Thus, the curves corresponding to that 
value of a would indicate those combinations of k s and ts for which 
the airplane is more stable with a given autopilot than without it. This 
value of a may be determined, of course, by solving for the complex 
roots of the characteristic stability equation when k s is set equal to 
zero. The family of curves obtained for this value of a forms the 
boundary between two different types of curves; that is, the curves for 
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dampings greater than the airplane damping are of a different shape than 
those for dampings less than this value. This change in the damping 
curves is discussed in detail in a subsequent paragraph. 


For purposes of illustrating the present method, the airplane is 
assumed to have freedom only in yaw. Several damping curves were calcu- 
lated for a typical present-day airplane equipped with an autopilot, 
sensitive to yawing acceleration for both the three-degree-of -freedom 
and the one-degree -of -freedom cases. The results of the calculations, 
although not always in quantitative agreement, indicated that the nature 
of the damping curves obtained for both cases was essentially the same 
and therefore subsequent discussion is based on the results obtained when 
the system was assumed to have freedom only in yaw. 


The mass and aerodynamic characteristics of the typical present-day 
airplane are presented in table I. The characteristic stability equa- 
tion for the three-degree-of -freedom case was presented previously 
(equation (4)), and the characteristic stability equation for the one- 
degree-of -freedom case is given by 

2p b K z 2 ! 2 - |c nr \ - = C n&i k s e '(7) 


The derivative C n ^ is assumed equal to -C n p. The equation used in the 
calculations, obtained by rearranging equation (7), is 


k s e 


■T S X 


2pb K z 2 ^ 2 - |C n r ^ 



( 8 ) 


As was pointed out in a previous section, the damping curve corre- 
sponding to the damping of the airplane oscillatory motion without the 
autopilot installed forms the boundary between two characteristically • 
different types of curves. In order to illustrate the types of curves ' 
which will be encountered, calculations were made, for the airplane 
damping and for dampings less than- and greater than this damping. The 
results are shown in figure 1. For the airplane described in table I, 
the airplane* damping is Tp /2 = 2.02 seconds. The values of t s and k s 
obtained from equation (8) have been converted from the nondimens ional 
time parameter s^ to time in seconds. Thus, the curves are plotted as 
a function of the dimensional parameters k and t in seconds. Each' 
point on these curves corresponds to a different frequency o) s , which 
varies from infinity at the point ( T=0,k=0. 0628) to zero as t .and k 
approach infinity. It is interesting to note that all the damping 
curves approach the point (0,0.0628) as cjd s — > 00 . For this point, the 
roots of the characteristic stability equation do not indicate any high- 
frequency oscillation; however, for any finite value of t, the roots of 
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this equation will indicate a high-frequency oscillation in addition to 
the usual lateral oscillation. For a given value of k, the frequency 
of this additional oscillation decreases with increases in r. Also, 
all the damping curves approach a minimum value of k for a value of o> s 
corresponding to a frequency close to the natural frequency of the air- 
plane. In general, frequencies greater than and less than the airplane 
natural frequency are separated by the point of minimum k on each of 
the damping curves. For this airplane the natural frequency is approxi- 
mately 5 radians per second. 

The curves presented in figure l(a) correspond to = °°. The 

complete boundary for neutral stability is made up of the curves for 
m = 0, 1, 2,. . . Unless m is allowed to take on integral values 
it is impossible to define correctly the stable region of the k,T plane. 
For example, consider the point (k=0.035,T=l. 6 ) in figure 1(a). The 
point, according to the curve for m = 1, represents a stable condition, 
but since it falls on the unstable side of the m = 2 curve, the pre- 
sence of unstable oscillatory modes in addition to the stable oscillation 
is indicated and thus this combination of k and t would result in an . 
unstable condition. The true boundary of neutral stability is shown as 
the solid curve of figure 1(a). Only the curves for m = 1, 2, and 3 
are shown' The curve for m = 0 would fall in the region corresponding 
to negative values of r and therefore was omitted. From this figure 
it is apparent that for some values of k there exist several ranges 
of r for which the system will be stable. The width of these ranges 
diminishes as m increases and eventually disappears completely. For 
this special case of T 1/2 = 00 and a given value of k, the ranges 
of t necessary for stability can be calculated analytically by the 
methods of reference 3, since the airplane was originally stable. The 
results obtained by using this method agree exactly with those of 
figure 1(a). The analysis presented in this reference was for a one- 
degree-of --freedom system and in appendix A it has been modified to apply 
to an airplane -autopilot system with three degrees of freedom. It should 
be emphasized that the method of reference 3 is applicable only to the 
special case for T 1/2 = 00 an ^ an originally stable airplane. On the 

other hand the method presented in this paper is equally applicable to 
originally unstable airplanes, in which cases the Tjy 2 = 00 curves will 

have the form of the curves -for dampings greater than the original air- 
plane damping. 

For T ^/2 = 3*50 seconds, the resulting curves for m = 0, 1, 2, 

and 3 are presented in figure l(b). The general shape of the curves is 
the same as for those of figure l(a), and the same -general conclusions 
apply. O nl y that part of the curve for m = 0 which is in the region 
of positive time lags is shown in the figure. The. curves for the air- 
plane damping, Tp/ 2 •= 2.02 seconds, are shown in figure l(c). For the 
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combinations of k and t within the triangular-like regions, which 
are formed by the curves for m, m + 1, and k = 0, the oscillatory 
motion of the an-piane is more stable with the autopilot than without 

e damping curves for T p/ 2 - 1.^0 seconds and TWo = 0.70 second 
are presented in figures 1(d) and 1(e), respectively, for m = 1- 2 
and 3 . These curves are quite different from the ones presented for 
damping less than the airplane damping. The combinations- of k and t 
for which the airplane would be more stable than is indicated by the - 
curves for T ± / 2 = 1.1*0 seconds or Ti/ 2 = 0.70 second, for exLple, 
are enclosed in island-like regions in the stable part of the k,r plane 

t^ llar l e fl° n exists for each va lue of m, but as m is increased 
the size of the region diminishes. In addition, as TWp is decreased 

these looped regions become increasingly smaller for each value of m 
As t 1/2 is reduced beyond some finite value these loops will cease to 
exist. The value of T x / 2 for which the m = 1 loop disappears 

represents the maximum damping obtainable for this airplane by the use 
of a specific autopilot. A method of obtaining this maximum damping . 
for each value of m is presented in appendix B. It can be seen from 
figure 1(e) that no loops exist for m > 1. This amount of damping, 
therefore, can only be obtained for all the oscillatory modes within 
the m - 1 loop. This conclusion could have been reached by carrVinc 
out the analysis of appendix B.. As a matter of interest, for this parti- 
cular autopilot -airplane system, the m = 1 loop ceases to exist for 
l l/2 < °*35 second. 


The regions in the k , t plane, which correspond to the dampings 
shown in figure 1, have been combined and are presented in figure 2. 

This figure gives an over-all picture of the effects of the specific 
autopilot on the lateral stability of the airplane. In order to verify 
the results indicated by this figure, airplane motions were calculated 
for the combinations of k and t identified as points A, B, C, and D. 
Step-by-step calculations were made by assuming freedom only in yaw, and 
the airplane motions in yaw subsequent to a 5 ° displacement in yaw for 
these combinations of k and t are presented in figure 3 . Point A 
is located within the m = 1 loop corresponding to T-^/g < l.ko seconds 

and the motion of figure 3(a) checks this prediction. The motion for 
point B (fig. 3(b)) is very lightly damped as is indicated by its near- 
ness to the zero-damping curve. Point C represents a combination of k 
and t for which the airplane is less stable with this autopilot than 
without it. Also, as seen from figure 2, the motion for this combina- 
tion of k and t should indicate two oscillatory modes of different 
frequency with approximately the same amount of damping. The motions- 
of figure 3 (c) support this prediction since more than one oscillatory 
mode is indicated. The combination of k and t for point D is located 
within the m = 2 loop, T 1/2 ^ 1.46 seconds. The motion for this combi- 
nation (fig. 3 (d)) is approximately as heavily damped as that of 
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figure 3(a), but because of the larger value of time lag this high 
damping is not realized as quickly. It is of significance, however, that 
this amount of damping can be obtained for the combinations of k and 
t within the m = 2, 3, . . . loops as well as for those within the 
m = 1 loop. In actual practice the phase lags encountered in conventional 
autopilots which can be represented by second-order systems generally do 
not exceed 180° at any frequency and the assumption of constant-time-lag 
characteristics usually is valid only within the frequency range where 
the autopilot phase lag is less than approximately 90°, The time lags 
common to the m = 2, 3, > • • loops imply phase shifts greater than 
2 it and thus are beyond the range where the assumption of constant time 
lag is usually valid. Two conceivable cases may be mentioned for which 
it might be desirable to develop devices with such large time lags. 

One is the case of a high-natural-frequency system (not necessarily an 
aerodynamic system) , for which' it might be impractical to obtain the 
small lags necessary for the stability of the system with automatic 
stabilization in the m = 1 region. Another case is that of an aero- 
dynamic system which is to respond to automatic control in addition to 
being automatically stabilized, in which case it might be desirable to 
delay the operation of the stabilization device. in order to reduce its 
effect on the response time of the automatic control. 

It should be reemphasized that the preceding analysis was made for 
Rfi airplane equipped with an autopilot having constant time lag. Damping 
curves such as those presented in figure 2 therefore apply, in a strict 
sense, only to autopilots which have these assumed characteristics. If, 
however, a frequency-response analysis of a given autopilot indicates 
that the assumption of a constant time lag is valid for only a limited 
range of frequencies, the preceding methods can still be used to good 
advantage. Damping curves such as those of figure 2 can be constructed 
and the values of k and t, which will result in a given amount of 
damping of the oscillatory modes of motion, can be determined since for 
each point in the k,r plane there will be certain complex solutions of 
the airplane -autopilot characteristic stability equation (equation (7)). 
If the frequencies of the oscillatory modes characteristic of a parti- 
cular point in the k,T plane are such that the autopilot can be made to 
exhibit constant-lag characteristics up to and beyond the highest 
frequency involved, then the damping of the oscillatory motions for the 
airplane -autopilot system will be essentially the same as predicted by 
the constant -time -lag analysis. The validity of this statement can be 
seen by examination of equation' (8) . The left-hand side of this equation 
is the autopilot frequency-response function & r /I>b Likewise the 
right-hand side of the equation is the airplane frequency-response func- 
tion br/Lb 2 ^- The necessar y 811(1 sufficient condition that X = a + ioig 

be a solution of the characteristic equation (equation ( 7 )) is that, 
when a + ico s is substituted for X in equation (8) , the complex number 
obtained from the autopilot frequency-response function be 
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identical with that obtained from the airplane frequency-response func- 
tion Sp/D^ 2 ^. The frequencies characteristic of any point in the k,T 

plane are readily identifiable and, if the frequency-response function 
Sr/Dt 2 * of an actual autopilot can be closely approximated by the ex- 

c — T X 

pression — r = k 8 e 8 , up to and beyond those frequencies 9 then the 
r>b 2 \)f 

complex roots of the airplane-autopilot characteristic equation must 
be essentially the same as those predicted by a constant-time-lag 
analysis. 


CONCLUDING REMARKS 


A method is presented for determining the control gearing and time 
lag necessary for a specified damping of an aircraft equipped with a 
constant-time-lag autopilot. The method is applied to a typical present- 
day airplane equipped with an autopilot which applies rudder control 
proportional to the yawing angular acceleration. The results calc ula ted 
for an airplane-autopilot system by the method described are in excellent 
agreement with the airplane motions calculated by a step-by-step procedure. 

The investigation shows that, for some values of control gearing, 
there exists more than one range of time lag for which the' assumed 
aircraft-autopilot system will have a specified amount of damp ing. The 
width of these ranges diminishes with increased time lags, and beyond 
some value of lag, these ranges cease to exist. 

It should be emphasized that the analysis presented in this paper 
is applicable, in a strict sense, only to systems which exhibit constant- 
time-lag characteristics; however, under certain conditions discussed 
in the paper, a constant-time -lag type of analysis is valid for systems 
which have constant time lags for only a limited range of frequencies. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , September 2 6 , 1950 
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APPENDIX A 

EXTENSION OF ANSOFF'S STABILITY CRITERION TO THE LATERAL 
MOTIONS OF AN AIRPLANE EQUIPPED WITH YAW AUTOPILOTS 


Reference 3 analyzes the effects of constant time lag in the feed- 
hack of a one-degree-of -freedom oscillatory system on the stability of 
the oscillation. The analysis shows that, in general, for a given 
gearing, ranges of time lag will exist for which the system remains 
stable. These ranges are separated from each other by unstable ranges. 
The expression for these stable time-lag ranges is obtained in closed 
form in terms of the other parameters of the system. It should be 
emphasized that the analysis is valid only for a system originally 
stable without feedback. 


This appendix extends the method of reference 3 to the three degrees 
of freedom of lateral motion of an airplane with a yaw autopilot. The 
argument and analysis presented in reference 3 are repeated herein only 
to the extent necessary to present the modifications of the method and 

to emphasize or clarify points. 

* 


The autopilot is assumed to provide a yawipg moment proportional 
to the yaw displacement or to the nth time derivative of the yaw, with 
the moment lagging the disturbance by a constant time t. The resulting 
nondime ns ional equations of motion are given by equations (l) and (2), 
and the characteristic stability equation of the system is equation (4). 
Note that the first five terms of equation (4) are the characteristic 
function of the airplane without the autopilot installed. If, therefore, 
equation (!<■) is divided through by this characteristic function of the 
airplane alone, the characteristic equation, of the total system is 
obtained in the form 


F(X) = 1 + 


(A’A.3 + BU 2 + CH D’) X n k 6 e" Tf?A ' 
AX 5 + BX 1 * + CX 3 + DX 2 + EX 


5 1 + W(X) =0 (Al) 


The theory may be given very briefly as follows: The roots of 

equation (4) are also the roots of equation (Al) and are the character- 
istic roots of the system. In order for the system to be stable, their 
real parts must all be negative. In the complex X-plane a closed curve 
consisting of the imaginary axis and the infinite semicircle to its 
right will enclose all values of X with positive real parts. Then by 
a theorem in the theory of functions of a complex variable (see refer- 
ence 6), it can be shown that the transformation of this closed curve 
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into the complex F-plane will circle the origin in the F-plane once in 
a positive direction for each root of F( X) =0 enclosed by the curve 
in the X-plane and once in the negative direction for each pole. It can 
be seen from equation (Al) that the origin of the F-plane corresponds 
to the point (-l,0i) in the W-plane, and, in fact, each curve in the 
F-plane is simply displaced unit distance to the left in the W-plane. 

If the airplane were originally stable, F( X) would have no poles in 
the right half of the X-plane; therefore, the condition that the whole 
system have no roots in this right half of the plane (unstable roots) 
is simply that W(X) . should not enclose the (-l,0i) point in the 
W-plane. 


One condition for the preceding theorem to apply is that there be 
no zeros or poles of F(X) on the boundary curve in the X-plane. It 
will be shown subsequently. that for n > 2 stability is impossible; 
thus, the only autopilots of interest are those where n = 0, 1, and 2, 
that is, displacement, velocity, or acceleration autopilots. Now the 
zero root of the characteristic equation (4) for n = 1 and 2 is 
removed by the division which gives equation (Al); therefore, no diffi- 
culty occurs for these two cases. For n = 0, however, equation (Al) 
gives a pole for F( X) at X = 0, which is on the boundary curve; 
therefore, for n = 0 dividing through by X is not convenient. Thus, 
for n = 0 


F(X) ■= X + W(X) = X + 



B'X 2 + C'X + D , )k s e~ TsX 
+ BX 3 + CX 2 + DX + E . 


= 0 


(A2) 


Finally, it is shown in reference 3 that, because of the negative 
exponential in equations (Al) and (A2), the transform in the W-plane of 
the boundary curve of the right half .of the X-plane degenerates to simply 
the transform of the imaginary axis. Thus, this curve can be obtained 
by placing im s for X in equations (Al) and (A2) and letting oo s 
vary from -00 to 00 . Equations (Al) and (A2) become 

F(io> s ) ■ 1 + W(itDg) = 0 . (A3) 

for acceleration and velocity autopilots, and 

F(ici)g) = ici) s + w(ia) s ) = 0 (a4) 

for displacement autopilots. 

For displacement autopilots, the critical point in the W-plane 
corresponding to the origin in the F-plane becomes the point -iu) s , a 
varying point; that is, W(a)g) must not enclose the point (0,-ia) s ). 
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Acceleration Autopilots 


The expression for W as a function of cDg may be obtained for 
the acceleration autopilot by use of equations (Al) and (A3). For the 
acceleration autopilot, letting a=A'k s , b=B'k s , c=C'k s , and 
d = D'k s gives 


W(ais) 


_( aCDs^ - CCDs 2 , 

+ 1 

^CDs " 

bco s 3) 

e - iT s a3 S 

(aod s ^ - Ca > s 2 + e) + i| 

(dco s - 

Bcd s 3) 


or 


W(u> s ) = 


(au) s ^ - cco s 2 ) 2 + (dco s - bcDg 3 ) 2 
^Ao> s ^ - Ccd s 2 + E^ 2 + ^Dcd s - Bco s 3^ 2 


1/2 


where 


0 = 0 2^ — Or 


e H-msr e ♦ 0) (a5) 


(A6) 


and 


6m s - boig 3 d - boig 2 

tan th = — _ 


axog' 1 ' - co3s^ aco s 5 - c <n s 


(A7a) 


tan = 


Da)s - Ba>s^ 
Ao>s^ - C (% 2 + E 


(A7b) 


The angles 6 ^ and $2 may be obtained from these equations for any 
given value of o) s since the algebraic sign of the numerator of tan 6 
is the algebraic sign of sin 6 and the algebraic sign of the denomi- 
nator is that of cos 0. These signs determine the proper quadrant 
for 0 . 

The W(a) s ) curve cannot encircle the point (-l,Oi) in the W-plane 
unless |w| >1. Thus, the ranges of cd s must be found where |w] > 1 
by setting |w| =1 and solving for (%. In the following discussion 
only the part of the curve for positive 0)3 will be considered inasmuch 
as the curve for negative u> s is simply the mirror image of the curve 
for positive o) s in the real axis. 

An important result can be obtained by considering what happens 
to W(cjd s ) as o) s 00 . The phase angle of W is seen to spiral around 
the origin whereas 
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lim IW 

(JDg — ^ 00 


|a| 



k s C n 5r 

|a| 

2n b (K x 2K z 2 - K XZ ) 

'V* 

2 ^b K z 2 


(a8) 


Thus, W(cd s ) spirals about a circle of this radius, and it can be seen 
that stability is impossible if the limit is greater than 1 , since the 
point (-l,Oi) would certainly be enclosed in this case. Therefore, 

stability is impossible for n > 2 since in this case lim Iwl = 

- co s — > 00 

In the cases of interest the limit of equation (A 8 ) is less than 1. 
The next step is to find the critical values of cn s by setting |w | = 1. 
This condition gives a quartic in cd s 2 . If x is substituted for <n s 2 , 
the critical roots are the positive real roots of the resulting equation 

(a 2 - A 2 -)x k + (b 2 - B 2 + 2AC - 2ac)x3 + (c 2 - C 2 + 2BD - 2 bd - 2 Ae)x 2 + 

(d 2 - D 2 + 2CE )x - E 2 = 0 (A9) 

When there are no positive real roots of equation (A9), the system 
is always stable. Moreover, for |^| < 1, there must always be an even 

number of such roots (double roots being counted twice). For the one- 
degree-of -freedom case considered in reference 3 this equation was a 
quadratic; consequently, there was only one pair of values (o> s ^,u> S 2 ) for 

which |w | =1, and between these values |w| >1. This condition corre- 
sponds to the fact that the frequency-response curve of the system has 
only one maximum. In the three -degree-of -freedom case there are four 
possible positive real roots, a condition which corresponds to the fact 
that the frequency-response curve may have two maximums (see fig. 4 (a)). 
Since either of these maximums, in case four critical roots exist, may 
enclose the point (-l,Oi), each pair of successive roots must be treated 
in the same way as the single pair was treated in reference 3 . In most 
cases equation (A9) will have only two positive real roots. A necessary 
condition that there be four positive roots is that equation (A9) have 
four sign changes in its coefficients. The necessary and sufficient 
conditions for existence of four real roots can be found in most algebra 
books ( for example , reference 7 ) • 

Although equation (A9) will usually have only one pair of real 
positive roots, for the sake of generality the possibility of two pairs 
must be considered. Corresponding to the two pairs of successive roots 
^cu Sl and a> S2 ^ and (o) Sl * and ), where u> S;L ' > afe ' > cn Sl > cd S2 , 

two pairs of values (h, * 2 ) and j$ 2 * ) can be obtained from equa- 

tions (A 6 ) and (A7). Then, there are two conditions to determine the 
stable ranges of time lag: 
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(2m-l)it + 0p (2m + l)n + 0T 

< t < 


0 ), 


s 2 


CD, 


'si 


(AlOa) 


(2m - 1 ) irt + 02 * 



< T < 


(2m + l)n + 0 1 ‘ 



(AlOb) 


For the system to be stable, both of these must be satisfied; that is, 
the stable time-lag ranges are those. where the ranges in equation (AlOa) 
and equation (AlOb) overlap. Here m = 0, 1, 2, . . . . The values 
of, m for which the stable, ranges exist are given by the condition 


or 


m < il t ^l) a)g 2 + ( n ~ 02 

2jt^cd S;l - cd S2 ^ 


(Alla) 


(* + 0i') oisg' + (* - 02 , ) c °si* 

2^cd S i ' - cd S2 ») 


(Allb) 


Of course the smaller of these two bounding values is actually the 
largest m for which stable ranges exist. This value can be used in 
the corresponding equation (AlOa) or (AlOb) to find the maximum value 
of t for stability. 


Velocity Autopilots 


The procedure for velocity autopilots is the. same as that for 
acceleration autopilots, the only difference being that the magnitude 
of W(o) s ) is divided by a) s and the phase angle decreased by it/2. 
The change in phase angle appears in 6±. f which is now given by 


tan 0]_ = 


co) s - aco s 3 
d - ba» s 2 


CA12) 


Setting | W | equal to 1 now yields, 
A 2 x^ + (b 2 - a 2 - 2 Ac)x3 + (c 2 

(d 2 - c 2 + 2b d - 2 Ce)x + E 2 


for u> s 2 = x, 

- b 2 + 2ac + 2AE - 2 Bd)x 2 

d 2 = 0 


+ 


(A13) 
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It might be noted that for this case lim Iwl = 0 

^s-* 00 ' 

similar to that of equation (A8) arises. However, 
does not vanish as a> s ->0: 


, so that no condition 
in this case, |w| 


U) c 


limJWl = 


d 

E 


2k s c Z)3 C n5 r 

Cn r Cz p -Cz r Cn p + (CZpCnp - C np Cz p )tan 7 


(Al4) 


Since |E| is generally very small, this quantity will generally be 
greater than 1. In this case, equation (A13) will have one or three 
critical roots (usually one), as may be seen from figures 4(b) and 4(c). 

In the case |j|| <1 an even number of critical roots will exist; and 

this case is treated exactly as in the case of the acceleration autopilot 
by using equations (A10) and the associated discussion, with equa- 
tion (A12) replacing equation (A7a). However, in the more general case 
of |j|| > 1, a slightly different treatment is needed. 


For the usual case of an odd number of critical values of co s ^ the 
lowest value is treated separately. Call this value oo SC) . Then, |w| >1 
when o> s < u) S q, as can be seen in figure 4(b). Since the phase angle 
( -o) s T + 0) is a monotonically decreasing function of and vanishes 

in this case for c% =0, .the first condition for stability becomes 

-^so 7 " + 00 > where 0o is obtained when o) S q is used in equa- 
tions (A12), (A7b), and (a 6). Thus, the stability condition for the 
first maximum of |w| is, the stringent one that 


t < 


* + 0Q 

“s 0 


(A15) 


In the case of one root, this is the only condition. For three roots an 
additional condition must be satisfied since the second maximum may - 
enclose the point (-l,0i). In this case the remaining two roots, with 
^Sq < ^82 < ^Sjj are paired just as in the acceleration autopilot, except 

that equation (A12) of course replaces (A7a), and equation (A10a) applies. 
Note that both conditions (equations (A10a) and (A15)) must be satisfied 
for stability, as described in the discussion of the acceleration autopilot. 


Displacement Autopilots 

In comparing equation (A2) with equation (Al) for n = 1, the W(a)g) 
in equation (a 4) for displacement autopilots can be seen to be exactly the 
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same as the W(to s ) in equation (A3) for velocity autopilots. The 
difference between the two cases arises from the change in the critical 
point of the W-plane, as pointed out in the paragraph following equa- 
tion ( a 4) . Only positive values of co s need be considered, and for 
these values the critical point is (0,-io>). Since this point cannot be 
enclosed except when |w| > cd s , the critical values of CDg 2 are given 
by the positive real roots of |w | 2 = cd s 2 , and the critical phase angles 
are now -^2m + rather than -(2m - l)n. 

The equation determining the critical values for x = <% 2 is 

A 2 x 5 + (B 2 - 2A0)x^ + (c 2 - a 2 + 2 AE - 2Bd)x3 + 

(d 2 - b^ + 2ac - 2CE)x 2 + (e 2 - - c 2 + 2bd)x - d 2 = 0 (Al6) 


This equation gives the intersections of the function |w| with the 
straight line with unit slope, as shown in figures 4(d) and 4(e). From 
these figures the equation is seen always to have an odd number of real 
positive roots (double roots being counted twice); and, in particular, 
figure 4(e) shows how the maximum of five critical values might occur. 


The method of procedure is the same as that given for the velocity 
autopilots with an odd number of roots. Equations (A12) , (>A7b) , and 
(A 6) are used to get 0(.o> s ). The first critical frequency a) So gives as 

the stability condition in this case that -<^so T + 0o > “o' ' Thus, the 


first stability condition for this case is apparently the very stringent 


one that 


§+0o 

^0 


(A17) 


For three critical values, the next two are paired- as before and 
give the additional stability condition for o) S q < a>s 2 < 


( 2 m - + 02 ( 2m + l)* + 01 

<% 2 “tei 


(A18) 


which replaces equation (AlOa).- Equations (A17) and (Al8) must both be 
satisfied for stability. 


Finally, for five critical values, an almost impossible case, the 
final pair of values ^s 2 ' ) gives a third necessary condition when 

used in equation (Al8), as described in the discussion of equation (AlOb). 
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APPENDIX B 

DETERMINATION OF MAXIMUM DAMPING OBTAINABLE BY USE 
OF A SPECIFIC AUTOPILOT 


The conditions which must he satisfied if the damping curves are to 
have loops are that ( T s) m have both a maximum and a minimum value for 
some values of u)g, say (o)g)^ and (<%)£, and that (kg) m be a mini- 
mum for some value of co 6 which is between ( 0 %)-,^ and (co s ) 2 * (See 

fig. 1(d).) The loops will break down when — ks ^ an d d ( T vanish 

dcn s do) s ‘ 

for the same value of. (o s . The value of a for which these conditions 
are satisfied is the maximum damping of the* system obtainable by use of 
a specific autopilot. A Curve representing a damping greater than this 
maximum in the k s /r s plane will have no practical significance, since 
for every point on this curve some oscillatory mode with less damping 
than the maximum will always exist. 


Now, (k s ) m and 


( T s) m are defined by the expressions 


u 6 ) m - 


( T s) m = 


2ran - Qr 


O)o 


(Bl) 


where R 2 , (Tg) m , and 6 2 are functions of cd s and were discussed in 

the section entitled "Analysis.” Differentiation of equations (Bl) 
results in the expressions 


d ( k s) m 

da)g 


( T s) r 


Ro e' — Rt s ) a] + e (Ts)ma ^ 

dm L J 


( s 


dlDc 


d ( T s) m 

do)g 


d0o 

“ 2rtm + 0 2 
o) s 2 


Setting these expressions equal to zero and considering m = 1 gives 
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dk 


— - = hR 2 — + — 

da) s diDg <i% 


= 0 


cLt, 

dio c 


s _ 


d 6 


(B2) 


= 2 * - 0 2 + <% eST - 0 


dr 


Since =0 is a necessary condition, ' 

cLd, 


dR 2 


must also "be equal to 


zero if < ^ ks is to vanish. • The expressions for R 2 and 0 2 are both 
<^s / \ 

functions of u>g, with the quantity a as a parameter. Equations (B2) 
should therefore he solved simultaneously in order to determine the 
maximum obtainable damping since 


T , - -0-69.3 * 

T l/2 " a V 


A convenient method for determining the maximum obtainable damping is 
as follows: For several values of a. for which the loops are known to 

be small, determine the values of tog for which = 0. Substitute 

UL ° S cLt s 

these combinations of a and u > 8 into the expression 3 -— and evaluate 

cLt s 

the function. A plot may then be made of against a as abscissa, 

and the value of a for which ^ = 0 is the required solution. - As • 

aia B dRo 

a check, this value of a should be substituted into — = 0 and the 

s 

corresponding value of cDg determined. This combination of a and oig, 
if a solution, will then satisfy ^ = 0 . 

The same analysis is applicable to the loops corresponding to 

m = 2 ,3 Equations (Bl) show that R 2 and 0 2 are independent 

dR^ 

of m; hence , 
for 


do, 


is also independent of m. The general expression 


d( T s)m _ q ( see equations (B2) ) 


dnig 


d 6 


- 0 2 + o) s ^ 


= 0 


Rearranging the equation gives 


d 0 o 

2ran = 0 2 ' dHT 


(B3) 
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Thus, if the combinations of a and which satisfy 


<3R 2 


are deter- 


mined for a sufficiently wide range of a and substituted into the 
right-hand side of equation (B3), the value of a for which loops will 


cease to exist for any value of m can be determined. For example, the 


m = 1 loops will break down when 0 2 “ ^s = the m = 2 loops, 

when this expression equals kit; and so on. In general, the value of a 
for which the loops break down decreases as m increases. 


The preceding analysis was applied to the m = 2 and m = 3 case 
for T 1/2 = 0.70 second since these loops were known to have broken 
down for this damping (see fig. 1(e)). As was expected, the analysis 
verified the fact that the loops did not exist for these values of m 
for T -jy 2 = 0.70 second. 
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TABLE I 


MASS CHARACTERISTICS AND STABILITY DERIVATIVES OF A 


W/S, lb/ sq ft 
S, sq ft . . 
b, ft . . . . 
p, slugs/cu ft 
V, ft/sec . . 
7, deg . . '. 

CL 

Pb 

..... 

Kz? 


K xz • 
Tb de s 


c l ’ 

per 

radian 

<V 

per 

radian 

C n p ^ 

per 

radian 

^nr-’ 

per 

radian 

C Y p , 

per 

radian 

CY r ; 

per 

radian 

°Yp. 

per 

radian 

c np> 

per 

radian 


per 

radian 


C n g > P er radian 


TYPICAL PRESENT-DAY AIRPLANE 


. . 65 

130 

28 

O.OOO89 
. 797 

. . 0 

. 0.23 

. 80.7 
0.00967 

0.0513 

-0.00145 

. - 2.0 
. -o.4o 

. 0.08 
-0.016 

. -0.4o 

. . 0 

. . 0 

. - 1.0 

. 0.25 
. -0.13 
-0.163 
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Acceleration autopilot 



Velocity autopilot 



Angular frequency, w 3 


Displacement autopilot 

Figure 4.- Types of frequency response for |w| showing possible critical 
frequencies for acceleration, velocity, and displacement yaw autopilots. 
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